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Abstract 

Stochastic difference equations and a stochastic partial differential equation (SPDE) 
are simultaneously derived for the time-dependent neutron angular density in a general 
three-dimensional medium where the neutron angular density is a function of position, 
direction, energy, and time. Special cases of the equations are given such as trans- 
port in one-dimensional plane geometry with isotropic scattering and transport in a 
homogeneous medium. The stochastic equations are derived from basic principles, i.e., 
from the changes that occur in a small time interval. Stochastic difference equations of 
the neutron angular density are constructed, taking into account the inherent random- 
ness in scatters, absorptions, and source neutrons. As the time interval decreases, the 
stochastic difference equations lead to a system of ltd stochastic differential equations 
(SDEs). As the energy, direction, and position intervals decrease, an SPDE is de- 
rived for the neutron angular density. Comparisons between numerical solutions of the 
stochastic difference equations and independently formulated Monte Carlo calculations 
support the accuracy of the derivations. 



Short running title: Stochastic Neutron Transport 

Keywords: stochastic partial differential equation; neutron transport equation; Ito system; 
stochastic model; Boltzmann transport equation. 

Mathematics subject classification: AMS(MOS) 82D75, 60H15, 82C70, 60H10, 65C30. 



1 



1 Introduction 



In the present investigation, stochastic versions of the deterministic neutron transport equa- 
tion are derived. Specifically, stochastic difference equations and a stochastic partial differ- 
ential equation (SPDE) are simultaneously derived that account for the random effects of 
absorptions, scatters, and source particles and generalize the standard deterministic neutron 
transport equation. Numerical approximations of the SPDE, through solution of the system 
of stochastic difference equations, provide approximations to the randomly varying neutron 
densities and yield insight into the random behavior of neutron transport. The stochastic 
transport equations are most useful for problems involving low numbers of neutrons. As 
the coefficient of variation is often approximately inversely proportional to the square root 
of the population size, the deterministic and stochastic transport equations yield essentially 
the same results for high numbers of neutrons. 

There are alternate but apparently equivalent ways to derive a system of stochastic 
differential equations (SDEs) for a randomly varying dynamical problem. The first way 
involves deriving a master equation for the random process [HJ [26]. A master equation is 
a differential form of the Chapman-Kolmogorov equation involving transition probabilities 
and is a probability conservation equation for the probabilities of separate states. If the 
transition densities are expanded in a parameter that defines the size of the fluctuations or 
jumps, then a Fokker-Planck equation or forward Kolmogorov equation is obtained in the 
first few terms of the expansion. As the probability density of an SDE system satisfies a 
certain forward Kolmogorov equation, this procedure infers a particular SDE system that 
approximates the random dynamics of the problem. In this procedure, a system of stochastic 
difference equations is not derived as an intermediate step. A second way to derive a system 
of SDEs for a randomly varying problem is by studying the changes in the process for a short 
time interval which gives a discrete stochastic model. The discrete stochastic model infers a 
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system of stochastic difference equations which, in turn, leads to an appropriate SDE system. 
For example, consider a randomly varying problem where X = [Xi, X2, ■ ■ ■ , Xn] t is a random 
vector of iV components for the problem. Let AX be the change in the process for a small 
time interval At. The expectations jl(X,t) = E(AX)/At and V{X,t) = E(AX(AX) T ) / At 
are determined and a stochastic difference equation approximation for the problem has the 
form: 

X(t + At) = X(t) + fI{X(t),t)At + (V(X(t), t)) 1/2 y/At ff t (1) 

where ff t is a vector of length N of independent normally distributed random numbers with 
zero mean and unit variance. Finally, the SDE system that approximates the behavior of 
the randomly varying process has the form 

dX{t) = fl{X, t)dt + (V(X, t)) 1 ' 2 dW(t) (2) 

where W(t) is a vector of length N of independent Wiener processes. It can be shown 
that the probability density of solutions of the stochastic system ([I]) or ([2]) approximates 
the probability density of the original randomly varying process [1, 4, 6J. Assume now that 
there are M possible changes in the process with probabilities pj At for j = 1, 2, . . . , M for 
small At. Also, assume that the jth change alters the ith component by Xj^. Then, the 
elements of /2 are given by /ij = Ylj=iPj^j,i an d the elements oi N x N matrix V are given 
by Vij = ^2jLiPj^j,Aj,i f° r i,l = 1,2, . . . N. In addition, an equivalent SDE system to (T2]) is 

dX(t) = jl(X, t)dt + C(X, t) dW*(t) (3) 

where N x M matrix C has elements Cjj = \^(p X - 2 and W*(t) is a vector of length M of 
independent Wiener processes. Furthermore, equation ((3]) is also obtained from the master- 
equation approach under the same assumptions |14J . Thus, the two derivation procedures 
produce identical SDE systems for this general A^-component and M-change process. 

These two derivation procedures produce very reasonable stochastic equation models for 
a given phenomenon. For randomly varying systems where the dependent variable depends 
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on time and on secondary independent variables, a stochastic partial differential equation 
(SPDE) may be derived by replacing the Wiener processes in the SDE system with appro- 
priate Brownian sheets and letting the intervals in the remaining independent variables go 
to zero. The resulting equation is an SPDE model for the phenomenon. 

In this paper, stochastic difference equations and a stochastic partial differential equation 
are derived for the transport of neutrons in matter. In neutron transport, captures, scatters, 
fissions, and source neutrons occur randomly. As a result, the neutron angular density varies 
stochastically. The relative magnitude of the random behavior of the angular density is 
pronounced for low neutron densities, such as during reactor startup, but decreases as the 
neutron density increases. The standard or deterministic neutron transport equation (or 
Boltzmann neutron transport equation) describes the expected or probable neutron angular 
density with respect to position, direction, energy, and time [9J. Solutions of the deterministic 
neutron transport equation provide average values of the neutron angular density; actual 
realizations with time of the neutron angular densities, that include random effects from 
neutron interactions and sources, are not obtained. 

A stochastic partial differential equation is derived in this paper for neutron transport in 
a general three-dimensional absorbing and anisotropic-scattering medium where the neutron 
angular density depends on position, direction, energy, and time. In the present investigation, 
the medium is assumed to be constant with respect to material composition, i.e., zero power 
noise. Special random effects, for example, from randomly varying boundary conditions or 
from a medium that is randomly varying [281 129] are not studied in the present investiga- 
tion although generalizations of the SPDE to approximate such conditions may be possible. 
Using the derived stochastic neutron transport equation, sample paths (realizations) of the 
randomly varying neutron angular densities can be approximately computed. In addition, 
after computing many sample paths, moments of the neutron densities, for example, can be 
estimated. The stochastic neutron transport equation is derived from basic principles, i.e., 
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from the changes in the system that occur in a small time interval. The dynamical system 
is studied to determine the different independent random changes that occur. Appropriate 
terms are identified for these changes in developing a stochastic difference system where all 
independent variables are discrete. As the time interval goes to zero, a certain stochastic 
differential equation (SDE) system is inferred (e.g., [U IH El [161 [25]). Next, multidimensional 
Brownian sheets replace the Wiener processes. As the intervals in the remaining indepen- 
dent variables go to zero, the SDE system leads to an SPDE [2J [3]. It is illustrated how the 
stochastic transport equation can be solved computationally through numerical solution of 
a stochastic difference system. 

The neutron transport equation is of fundamental importance in nuclear reactor theory 
and shielding design P [121 El 122]. The stochastic nature of the neutron transport process 
has been of interest for many years. Classic studies of the stochastic theory of neutron 
transport are given in E]. In particular, let p n (R,tf;x,v,t) be the probability that a 
neutron with position x and velocity v at time t leads to n neutrons in region R of x , v space 
at time £/. In |8J, a non-linear integro-differential equation for the probability generating 
function is derived for p n (R,tf] x, v, t). The equations derived are interesting but complicated 
and difficult to apply. More recently, a master equation approach was used to estimate the 
temporal evolution of the number of neutrons in time- varying multiplying systems [19], [23] . 
This approach gives, for example, moments of the number of neutrons in the system with 
time. However, a stochastic difference system approximation of neutron transport is not 
determined and, as a result, sample paths of the randomly varying neutron densities with 
respect to energy, position, and direction are not estimated. 

In the next section, stochastic difference equations and a stochastic partial differential 
equation are derived for neutron transport in general three-dimensional xyz-geometry. The 
changes due to absorptions, fissions, and scatters, which occur randomly with probability 
proportional to the neutron angular density and to the time interval, are carefully consid- 
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ered in deriving the equations. In the third section, several special but useful cases of the 
stochastic neutron transport equation are described such as one-dimensional slab geometry 
with isotropic scattering. In the fourth section, stochastic difference equations are solved 
computationally for the randomly varying neutron densities and compared with Monte Carlo 
calculations. The Monte Carlo calculational procedures differ considerably from the numer- 
ical solution of the stochastic transport equations. In the Monte Carlo calculations, the 
dynamical system is checked at each small interval of time to take into account scatters, ab- 
sorptions, and movements for individual neutrons. Comparisons between the two different 
computational methods are in close agreement indicating that the stochastic difference equa- 
tions and the stochastic partial differential equation accurately model the random behavior 
of neutron transport. Although the stochastic equations cannot exactly model the random 
behavior of neutron transport since, for example, the number of neutrons in any region is 
not integer-valued in the stochastic difference system or in the SPDE model, the calculations 
indicate that the stochastic neutron transport model is accurate. In addition, the stochas- 
tic transport equations provide insight into the random dynamics of neutron transport and 
can be efficiently solved computationally using the stochastic difference equations. Several 
applications of the stochastic neutron transport equation are discussed in the fifth section 
before the investigation is summarized in the final section. 
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2 Derivation of Stochastic Neutron Transport Equa- 
tions 



The neutron transport equation in xyz-geometry can be written in the integro-differential 
form 0QJ2]: 

dN(x,y,z,u,<f>,E,t) dN ON dN 

dt = ~ V ^ ~ V ^ (4) 

+ Q(x, y, z, E, t) - va(x, y, z, E)N(x, y, z, fi, 0, E) 

/•Emax pi p^TT 

+ I / / v'a' f(x,y,z,n' ,<j)' ,E' ,[i,(j),E)N' dftdji'dE' 



J-UO 
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for (x,y,z) E ([0 

5 ^moi] X [0, ymax\ X [0,z max ]) C W, -1 < fi < 1, < < 2tt, < 
E < E max and t > where iV = N(x,y, z, fi, 0, E,t) is the expected neutron angu- 
lar density with respect to position (x,y,z), direction Q = (//,</>), energy E, and time 
t per unit volume per unit solid angle per unit energy. In (jlj, a' = a(x,y, z, E') and 
N' = N(x, y, z, fi', cf)' , E' , t). The notation used here is generally consistent with the notation 
used in [9J. In particular, v is the neutron speed, a is the total macroscopic cross section, 
a(x, y, z, E')f(x, y, z, fi' , (ft', E', fi, 0, E)A/j l A<f)AE is the probability of a neutron transfer from 
direction (fi', 0') and energy E' to solid angle A/iA0 about direction (fi, 0) with energy AE 
about energy E. Note that 

°'f(x, V, z, fi', 0', E' , fi, <f),E) = Y^ KMx, y, z, fi', 0', E', fi, 0, E) 

r 

where a' r = a r (x, y, z, E') and the sum includes the separate interactions r in which neutrons 
are produced such as elastic scattering or fission. The parameters fi x , fi y and fi z are the 

direction cosines for the x, y, and z axes, respectively. In particular, fi x = fi = cos 9, 

- - - dN dN dN 
u v = sinfcos -1 u) cos 0, and u z = sinfcos" 1 a) sin0. Thus, Q-VN = u x -^ — h u v -z — hMz^— • 

ox ay oz 

Furthermore, c(x, y, z, E) = j^ max j\ f(x, y, z, fi, 0, E, fi', 0', E') d(f)'dfi'dE' is the mean 
number of neutrons emerging per collision of neutrons of energy E at position (x,y,z). 
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Finally, Q(x,y, z, E) is the number of source neutrons per unit solid angle per unit volume 
per unit energy and is assumed to be isotropic. 

To simplify the derivation, it is useful to define several other quantities. Let a c be the 
capture cross section, i.e., the sum of all the cross sections involving a pure capture event 
such as those due to (n, 7), (n,p), or (n, a) collisions. Let a(x,y, z, E) = a(x,y, z, E) — 
a c (x,y,z,E) be the macroscopic cross section for all interactions other than pure capture 
interactions. Let / be defined by the expression 

a(x, y, z, E')f(x, y, z, //, <p\ E', fi, <fi, E) = a'f(x, y, z, //, 0', E r , //, 0, E) 

and define c(x, y, z, E) = J Q max j\ f(x, y, z, /i, 0, E, //, <//, E') dtp'd/i'dE' as the mean 
number of neutrons emerging per non-capture collision of neutrons of energy E at position 
(x,y,z). 

Equation Q is deterministic and random variations in the neutron angular density due to 
the inherent randomness in absorptions, scatters, and fissions cannot be accurately studied 
using this equation. To derive a stochastic partial differential equation generalization of (]1J), 
the changes which occur in the angular density for a small time interval are determined 
taking into account interactions and transport. A discrete stochastic model of the neutron 
angular density is then constructed which infers a system of stochastic difference equations. 
As the time interval decreases, the stochastic difference system leads to a system of Ito 
stochastic differential equations. As the intervals in position, direction, and energy decrease, 
a stochastic partial differential equation is derived for the neutron transport process. 

2.1 Several Properties of Brownian Sheets 

Before deriving these stochastic equations, it is useful to consider several properties of Brow- 
nian sheets [5j [101 E2]. A Brownian sheet on [0,5] x [0,5] is illustrated in Fig. [TJ The 
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Figure 1: A Brownian sheet on [0,5] x [0,5]. 



Brownian sheet W(x,t) satisfies: 

rt+At rx+Ax q2 W ( x > f) 

That is, the Brownian sheet is independent and normally distributed over rectangular regions. 
In addition, if Xj = jAx for j — 0, 1, . . . , J, where Ax = x max /J, then the Brownian sheet 
defines for j = 1, 2, ... J, the standard Wiener processes, Wj(t), where 

dWAt) 1 f** d 2 W(x',t) , . 
' ' — ' dx . 



dt v / Ax Jx--i dtdx 
Notice that if U = iAt for i — 0, 1, . . . , M, then 



r dw J (t , )=VAt Vi j 

Jti-x 



where rjij ~ A/"(0, 1) for each j = 1,2, ... J and i = 1, 2, . . . M. Also, standard Wiener 
processes can be defined using, for example, three-dimensional Brownian sheets letting 

/ x 1 f Xj f Vk d 3 W(x',y',t) , , , . , 



VAxAy } x ._ x J Vk _ x dtdy'dx 1 
where Wj t k(t) is a Wiener process for each j and k. However, notice that 

d 3 w(x, y ,t) = Um i r x+Ax r y+Ay d 3 w(x', y ',t) d , dx , 

dtdydx Ax,Ay^o Ax Ay J x J dtdy'dx' 



Finally, it is useful to note that W(x,t) ^ W 1 (x)W 2 (t) where W\{x) and W 2 (t) are 
independent Wiener processes. To see this, let A = [a,b] x [c,d\ and \A\ — (d — c)(b — a). 
Then, 

W(A)= ( [ ' - dxdy = W(b, d) - W{b, c) - W(a, d) + W(a, c) ~ MO, |A|). 

In particular, = rj xt where rj xt ~ A/"(0, 1). However, 

W h2 {A) = f f dW 1 (x)dW 2 (t) = (W^b) - W 1 (a)){W 2 (d) - W 2 {c)) = ^/\A\ Vx r] t 

J a J c 

where rj x ~ Af(0, 1) and ^ ~ jV(0, 1). So, for example, £((W(A)) 4 ) = 3|A| 2 whereas 
E{{W h2 {A)Y) = 9|^| 2 - 

2.2 Derivation of stochastic neutron transport equations 

To derive a stochastic neutron transport equation, the changes which occur in the neu- 
tron angular density for a small time interval At at time t p are considered where t p = 
(p — 1) At for p = 1,2, ... . To facilitate finding these changes, the variables position, di- 
rection, and energy are made discrete. Three-dimensional space is discretized into rect- 
angular parallelepipeds of length, width, and height Ax, Ay, Az, respectively. The di- 
rection variables p, and are discretized with intervals of size Ap = 2/L and A<p = 
2ir/M, and energy is discretized into intervals of size AE = E max /G. Furthermore, Xi = 
(i - I) Ax for i = 1, 2, ...,/ + 1, y 3 = (j - I) Ay for j = 1,2,. . . ,J + 1, z k = (k - 
l)Az for k = 1, 2, . . . , K + 1, p l = -1 + (/ - l)Ap for I = 1, 2, . . . , L + 1, <p m = (m - 
1)A0 form = 1, 2, . . . , M + 1, and E g = (g - 1)AE for g = 1, 2, . . . , G + 1. Now let 
n(xi,yj, z k , pi,<p m , Eg,t p ) = N(x i ,y j ,z k ,pi,(j) m ,E g ,t p )AxAyAzApA(j)AE be the number of 
neutrons in the parallelepiped moving in direction pi,(f> m with energy E g . There are several 
possible changes that can occur to n — n(xi,yj, z k , pi, (fi m , E g ,t p ) in the small time interval 
At. A capture or fission can occur, a neutron can enter or leave one of the six faces of the 
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parallelepiped, or a scatter can occur resulting in a loss or gain of one neutron. The possible 
changes An along with their probabilities are listed in Table 1 for a small time interval At. 
Notice that position changes occur deterministically, i.e., the neutron position is determined 
by the neutron velocity and, thus, the number of neutrons moving from one parallelepiped 
into an adjacent parallelepiped is calculated based on the fraction of neutrons crossing the 
parallelepiped boundary in time At. In Table 1, the probabilities p c ,Ptri,Ptr2, and pq are 
given by 

Pc = v g a c n(x i ,y j ,Zk,fJ,i,<l>m,Eg,tp)At 
Ptri = v g crfnAfiA(j)AEAt, 
Ptr2 = Vg'cx' f'n AfiA<f)AEAt, and 

p Q = Q(xi,yj,z k ,E g ,t p )AxAyAzAfiA(f)AEAt. 

For example, in Table 1, v g a c (xi, yj, Zk, E g )n(xi, yj, z k , Hi, <f> m , E g , t p )At is the probability that 
a neutron in the packet undergoes a capture. However, the two interaction terms, i.e., the 
term involving transfer into the packet and the term involving transfer out of the packet, need 
to be considered at each position for all the different directions and energies. Specifically, 

Ptri = v g af{xi, y^ z k , n h <f> m , E g , n v , 4> m >, E g >)nA/j,A(pAEAt 

is the probability in time At that l/c(xi,yj, Zk, E g ) neutrons are lost from direction u;,0 m 
at energy E g (from the packet) when one neutron emerges in direction /ii>,(f) m > with energy 
Eg/ for each value of /', ml, and g' and 

Ptr2 = v g/ a'f(xi, yj, z k , \i v , <p m ,, E g/ ,n h <p m , E g )n' AuA0A£"At 

is the probability that one neutron emerges in direction m at energy E g (into the packet) 
when l/c(xi,yj, Zk, E g i) neutrons are lost from direction fii>,(fi m > and energy E g i for each 
value of I', m' , and g' . Note, for convenience in Table 1, n — n(xi, yj, Zk, Hi, 4> m , E g , t p ), n' = 
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n(xi, i/j, Zk, Hi', 4> m ', Eg',t p ), a = a(xi, yj, z k , E g ), a' = a(xi, yj, z k , E g >), a c = a c (%i, Vj, z k , E g ), 

a = cr(xi,yj,z k ,E g ), f = f{x h yj,z k , Hi,<j) m ,Eg, HvAm' ,E g >), f = f{x h y h z k , Hi',4>m', E g >, Hi,4>m, E g ), 

and c = c(xi,yj, z k , E g ). Also, it is assumed that the neutron source, Q, is a Poisson process 

with the probability of adding one source neutron to the packet in a small time interval At 

equal to Q(xi, yj, z k , E g , t p )AxAyAzAfiA(f)AEAt. 

Table 1 defines a discrete stochastic model for the neutron transport system. Using these 
changes and probabilities and letting the time interval At approach zero, a system of Ito 
stochastic differential equations can be formulated for the dynamics of this random transport 
process. First, a deterministic equation for the expected number of neutrons at time t p + At 
can be derived using the results of the Table 1. This equation, for Hx, Hy, /i z > 0, is given by: 

n(xi, yj, z k , Hi, 4> m , E g , t p + At) = n(xi, y jy z k , Hi, <f) m , E g , t p ) (5) 
-v g <j c n(xi, y jy z k , Hi, 4> m , E g , t p )At 

+HxV g At(n(x i -i,y j , z k , Hi, 4> m , E g , t p ) - n(xi, y j} z k , hi, <f) m , E g , t p ))/Ax 
+fi y VgAt(n(xi, yj_ h z k , Hi, <fi m , E g , t p ) - n(xi, yj, z k , hi, 4>m, E g , t p ))/Ay 
+ k i z v g At(n(x i , yj, z k -\, Hi, <j) m , E g , t p ) - n(xi, yj, z k , Hi, 4> m , E g , t p ))/Az 
+Q(xi, yj, Zk, Eg, t p )AxAyAzAfiA<j)AEAt 

L M G j 

- ^2 - v a°H x i, Vj, z k, Hi, 4>m, Eg, flu, 4>m>, Eg^nAnA^AE At 

U=l m '=l g'=l 
L M G 

+ ^2 ^2 v g' a 'f( x i' Vh z k, Hi', <firn', E g ,, Hi, <fim, E g )n' AnA(j)AEAt 

l' = l m' = l g' = l 

where n is the expected number of neutrons in the packet. Letting n = NAxAyAzAnA<f)AE 
and allowing At to approach zero as well as Ax, Ay, Az, Ah, A<p, and AE, it is straight- 
forward to show that ([5]) yields the standard neutron transport equation (jlj). However, the 
changes and probabilities given in Table 1 can be used to derive an SDE model. Indeed, the 
discrete stochastic model and the SDE model will have approximately the same covariance 
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terms as well as mean terms for small At. 

The derivation procedure, described in the Introduction for obtaining equation (JHJ), is now 
applied using the changes and probabilities given in Table 1 to obtain the stochastic terms 
in the equations. Specifically, for the zth component (packet) of the system, the coefficient 
of the independent Wiener process corresponding to the jth change is equal to the product 
of \jj with the square root of the probability for the jth change, recalling that A^j is the 
amount that the jth change alters the ith component. The changes and probabilities given 
in Table 1 imply, for [i x ,{i y ,{i z > 0, that a very reasonable approximation to the discrete 
stochastic model satisfies the stochastic difference system: 

n(x h yj, z k , Hi, <f) m , E g , t + At) = n(x h y h z k , Hi, <f) m , E g , t) (6) 
+fi x v g (n(x i _ 1 ,y j , z k , hi, <p m , E g , t) - n{xi, y h z k , Hi, <t> m , E g , t)) At/ Ax 
+HyV g {n(xi, z k , Hh 4>m, E g , t) - n(xi, yj, z k , hi, 4> m , E g , t)) At/ Ay 
+HzV g {n{xi, y^ z k -i, Hh <fim, E g , t) - n{x h y d , z k , Hi, <j> m , E g , t)) At/ Az 
+Q{x u yj, z k , E g , t)AxAyAzAnA(f)AEAt - v g a c n(xi, yj, z k , Hi, <j> m , E g , t)At 

L M G j 

^ ^ ^-Vgaf(xi,yj,z k ,Hi,(f) m ,E g ,Hii,<j)rni,E g ,)nAHA(j)AEAt 

l'=l m'=l g'=l 
L M G 

+ ^ v 9' a 'f( x h Vj, z k , Hi', <Pm>, E g i,Hi, (frm, E g )n'AHA(f)AEAt 

l'=l m '=l g'=l 



+ \ Q(xi, yj, z k ,E g , t)AxAyAzAHA(pAEAt r] (Q) 



i,j,k,l,m,g 



(c) 

i,j,k,l,m,g 



-\Jv g a c n(xi, yj, z k , hi, 4> m , E g , t)At r), 

L m g 

"EEE -JvgafnAnA^AEAt rj^ 



c 

l'=l m'=l g'=l 
L M G 

+ EEE y/v^f'n'A^AEAt vSk, W ,i, m , g 

l'=l m '=l g'=l 

for i = 1,2,...,/, j = 1,2,..., J, k = 1,2,..., K, I = 1,2,..., L, m = 1,2,..., M, 
and g = 1,2,..., G, where vl%,i,m, 9 > V^mm and 4jl,i, m , g ,i',m', g ' are independent nor- 
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mally distributed numbers with zero mean and unit variance processes for each value of 
i,j,k,l,m,g,l',m',g'. In Equation ©, / = f{xi,yj,z k ,Hi,4>m,E g ,Hi',4>muE g ') and /' = 
f(xi,yj, Zk, Hi'-. 4> m ', E g ', Hi, <p m , Eg). For small At, the stochastic difference system (JSJ) has 
the same mean and mean square changes as the discrete stochastic model defined by Table 
1. 

Stochastic difference system fl6]) is an Euler-Maruyama approximation to a certain Ito 
SDE system (e.g., [U HJ E]) which has the form: 

dn(xi,yj,z k ,Hi,4>m,E g ,t) __ 



dt 

+HxV g (n(x i _ 1 ,y j , z k , hi, <frm, E g , t) - n(xi, y jt z k , Hi, 4>m, E g , t))/Ax 
+HyV g (n(xi, yj-i, z k , hi, 4>m, E g , t) - n(x i: y j} z k , hi, 4> m , E g , t))/Ay 
+HzV g (n(xi, y j} z k ~i, Hu 4>m, E g , t) - n(xi,yj, z k , hi, 4> m , E g , t))/Az 
+Q{xi, yj, z k ,E g , t)AxAyAzAnA(pAE - v g a c n(xi, y jt z k , Hi,<f>m, E g , t) 

L M G - 

^2-v g arf(x i ,y j ,z k ,Hi,4>m,Eg,Hi', ( Pm',E g ')nAHA(f)AE 

l'=l m '=l g'=l 
L M G 

+ v g' a 'f( x u Vj, Zk, Hi 1 , 0m', E g ,, hi, (frm, E g )n'AHA(j)AE 

l'=l m '=l g'=l 

dWgl, m /t) 
dt 



Q{xi, yj, z k , E g ,t)AxAyAzAnA(pAE 
-yj v g a c n{Xi,yj, z k , Hi,(p m , E g ,t) — 

L M G , , rrxAtr) 



dW { ' (t) 

aVV i,j,k,l,m,g,l'm',g'\ L ) 



EEE \\jv g afnAHA4>AE 

l'=l m '=l g'=l 



L M G dw (tr) (f) 

UV]/ i,j,k,l' \m' ',g' ,l,m,g\ L > 

dt 



+ Yl J2 s/vgxr'f'n'&n^E 

l'=l m' = l g'=l 



(7) 



for i = 1,2,...,/, j = 1,2,..., J, k = 1,2,. ..,K, I = 1,2,. . . , L, m = 1,2,. . . , M, 
and g = 1,2,..., G, where W t %^ g (t), W t %^ g (t), and \Y^, J>) are inde- 
pendent Wiener processes for each value of i, j,k,l,m, g,V ,m' , g' . In Equation ([7]), / = 
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f(x i ,yj,Zk,fJ>i,<l>rn,Eg,tJ>i',<l>m>,Eg') and f = f{x h yj, z k , Hi' , (f) m > , E g > , fj, h (f> m , E g ) . For small 
At, the stochastic system ([7]) has approximately the same mean and mean square changes 
as the discrete stochastic model defined by Table 1. 

Before the intervals in space, energy, and direction can be allowed to go to zero so 
that the SDE system will approach an SPDE, the Wiener processes need to be replaced 
with appropriate Brownian sheets. Introduced now are multidimensional Brownian sheets 
(x, y, z, //, 4>, E, t), W^(x, y, z, /i, <fi, E, t), and (x,y,z,n,</>,E,fj?,<f/,E',t). For ex- 
ample, W <yC \x,y, z, /i, 0, E,t) is an independent seven-dimensional Brownian sheet in vari- 
ables x, y, z, fi, 4>, E, t. The Wiener processes in (J7|) are now replaced by equivalent forms 
involving Brownian sheets after which the spatial, angular, and energy intervals will be 
allowed to approach zero. Specifically, in (J7J), let 

dW $l,l,m,g(t) 1 [ Xl+1 f Vj+1 f Zk+1 /" / ' m+1 



dt ^AxAyAzAfiA^AE J x 



E, 



9+ 



1 d 7 W^(x,y,z,HA,E,t) 



dE d(fi d/i dz dy dx 



dxdydzdfid(f)dEdt 
dW-\ (t) 1 r^ 1 f Zk +^ r^ 1 

^ ' * t,3,k,l,m,g\ u J _ 1 I I I I 



dt y/AxAyAzAiiAcpAE 

E a +i d 7 W^{x,y,z,n,<P,E,t) 



dxdydzd[id(f)dEdt 



dE dip dfi dz dy dx 



dW%i Vm ,„,r m „{t) l rw+i f". + i w +1 r<t> m > + i rE g/+1 



a 



dt VAS j Xi j Vj j Zk j H , J4, m , je, 

r<p m+ i rE g+1 d^w {tr) (x,y,z,u!,(j)',E',a,4),E,t) Irn 
' ' dE ■ ■ ■ dx 



m J(j>. 



dxdydzdfx'd(f)'dE'djjd(f)dEdt 

where AS = AxAyAzAfi'A<p'AE'AfiA<f)AE. These equivalent expressions are now substi- 
tuted into (J?D and n is replaced with NAxAyAzAfiA<pAE. Next, Ax, Ay, Az, A/i, A<p, and 
AE are allowed to approach zero. The result is a stochastic partial differential equation for 



15 



stochastic neutron transport: 




8 



+ / / / va'f(x,y,z,n',(f)',E',^(j),E)N'd(j)'diJ 1 ; dE' 



Jo J -i Jo 




dE'dQ 



dE'dVL' 



where, for convenience, N = N(x, y, z, fi, 0, E, t), N' = N(x, y, z, //, 0', E', t), a = a(x, y, z, E), 
a' = cr(x,y,z,E'), a(x,y,z,E) = a(x,y,z,E) - a c (x, y, z, E), c = c(x,y,z,E), [i x = /i, 
fiy = sin(cos _1 fi) cos0, fi z = sin(cos _1 /i) sin0, dVL' = dtp'dfj,', and the equation is valid 
for positive or negative values of the direction cosines fi x ,fi y ,fi z . Notice that N in (JS} is 
stochastic, i.e., iV is not the expected neutron angular density. That is, each solution of 
OH]) is one possible realization or sample path of the neutron angular density. There are an 
infinite number of these sample path solutions, as occurs in nature, and the average of these 
random solutions is equal to the expected neutron angular density N. Finally, notice that 
([8]) generalizes (TjJ. If the stochastic terms are set equal to zero, then ([8]) is identical to @. 
Of course, for the deterministic or the stochastic version of the neutron transport equation, 
the initial condition and boundary conditions must be specified. 

As illustrated in Section @] for two special cases of OH]), the stochastic transport equation 
(JH1) can be solved computationally by discretizing position, direction, and energy and then 
approximating the resulting system of Ito stochastic differential equations in time. In effect, 
SDE system (JTj) is computationally solved using an appropriate stochastic difference system 
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such as (J6j). 



3 Special Cases of the Stochastic Neutron Transport 
Equation 

In this section, for illustrative purposes, two special cases of (jHJ) are considered. First, a one- 
dimensional stochastic neutron transport equation with isotropic scattering is described. 
Second, a stochastic partial differential equation is given for neutrons interacting in a homo- 
geneous medium. 

Consider the parallelepiped region < x < x max , < y < y ma x, and < z < z max where 
the medium is uniform with respect to the spatial variables y and z. Assume that, except 
for the left face and right face defined by x = and x = x max , respectively, the neutrons 
are reflected back at the other faces, i.e., there are reflecting boundary conditions at the 
four faces except for the left and right faces. Equation (JHJ) is integrated over < y < y ma x, 
< z < z max , and angle (p. Then, the equation reduces to 

dN(x,u,E,t) dN ,_, . , , T . ^ . , n . 

L = ~ V ^1T + Q x > * - va ( x > E ) N ( X > ^ E > *) 9 

at ax 

/Emax /"l 
/ v'a(x, E')f(x, //, E', E)N(x, //, E', t) d^'dE 1 

m ^ w^\x,^E,t) — ^aV^^ 



1 Jo 



+ 



c(x,E) dxdfidEdfi'dE'dt 
7_ 1 7o v J dxdyfdE'dndEdt P 



where N(x, //, E, t) is the number of neutrons per unit length per unit angle per unit energy at 
position x with direction /x and energy E at time t, c(x, E) = j Q Emax j\ f(x, fi, E, //, E')d/j,'dE' 
is the mean number of neutrons emerging per non-capture collision of neutrons of energy E 
at position x, N' = N(x, //, E', t), and a' = a(x, E'). In Equation (jHJ), /' = f(x, //, E r , fi, E) 
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and / = f(x,fi,E,fi',E'). Equation fl9]) is a stochastic neutron transport equation for one- 
dimensional slab geometry. 

Furthermore, assuming a single energy, isotropic scattering, and only capture and scat- 
tering interactions, the above equation becomes: 

dN(x, /i, t) 



dt 



ON 1 f 1 

-va— h Q — va(x)N H — / va s (x)N(x, [i , t) d/J,' 

ox 2 J_ x 



(10) 



where A" = N(x,n,t), Q = Q(x,t), and cr s (x) is the scattering cross section at position x. 
Equation ({TO]) is a stochastic neutron transport equation for mono-energetic transport in 
one-dimensional plane geometry with isotropic scattering. 

Consider again the parallelepiped region < x < x max , < y < y ma x, and < z < z max 
where the medium is uniform with respect to all the spatial variables x, y, and z. Assume that 
the neutrons are reflected back at all six faces, i.e., there are reflecting boundary conditions 
at all the faces. Equation ([8]) is integrated over the volume < x < x max , < y < y ma x, 
and < z < z max , and over the angles 4> and \i. Then, equation ([H]) reduces to 

v'a(E')f(E',E)N(E',t) dE' - ^Jva c {E)N ^ ' ' 

- / -^ T y/va(E)f(E,E')N(E,t) \ \ V ' ' dE' 
J c(Ey ; v ' ; OEdE'dt 

where A" = N(E,t) is the number of neutrons per unit energy and Q(E,t) is equal to the 
number of source neutrons per unit energy per unit time. 

In a homogeneous medium, with only capture and scattering interactions, the above 
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equation becomes: 

^^ va{E)N + Q(E , t) + ^^^EA (12) 

v'a(E')f(E',E)N(E',t) dE' - a/ va c (E)N ^ ' ) 

Jva(E)f(E,E')N(E,t) \ | ' ' J dE' 
V v jj\ , ) K > j QEdE'dt 

where v'a(E')f(E',E) is the probability of a neutron scattering from energy E' to E per 
unit energy per unit time. Equation f fl2|) is a stochastic neutron transport equation for a 
homogeneous medium with only capture and scattering collisions. 

4 Comparison With Monte Carlo Calculations 

In this section, the stochastic difference equations derived in the previous sections for neutron 
transport are numerically solved and compared with independent Monte Carlo computations. 
Two cases are considered. First, a problem is studied involving mono-energetic neutron 
transport in a slab with captures and isotropic scatters. The stochastic transport equation 
in this case is given by equation ( 110]) . Second, a homogeneous medium is studied where 
the neutrons experience captures or scatters with energy changes. The stochastic transport 
equation for this problem is ( fl2l) . 

In the first problem, it is assumed that 1000 neutrons per second begin entering the left 
side of a slab of width one unit at time t = 0.0. The velocity of the neutrons is v — 0.1. The 
slab is homogeneous and the scattering and capture cross sections in the slab are assumed 
to be a s = 5.0 and o c = 0.10. Therefore, for this problem, the slab has width l/a c = 10 
absorption mean free paths. (As slab width increases, the leakage decreases but it is likely 
that the coefficient of variation in the leakage increases.) The neutrons isotropically enter 
the slab on the left side, x = 0, from time t = 0.0 until t = 50.0. After time t = 50.0, the 
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neutrons no longer enter the slab and the neutrons in the slab eventually are absorbed or 
escape. There is no neutron source for this problem. The time dependence of the exiting 
fluxes on the left and right sides of the slab are of interest in this problem. To study this 
problem computationally, equation f JTOj) needs to be solved numerically. 

To define a numerical method for equation (fTUl) . the interval [0, 1] in x is divided into / 
intervals [xj, Xj+i], for % — 1, 2, . . . , / where Xi — (i — l)Ax, and Ax = 1/1. In addition, the 
interval [—1, 1] in direction \i is divided into J equal intervals of width Afj, = 2/ J and time is 
discretized where t k = kAt, k = 0, 1, 2, ... . Considering ( flCfl) at position Xi and direction fij 
and using an explicit approximation in time t along with an upwind differencing approach 
suggests the numerical procedure: 

( HjVrii^j^At/Ax - fijVn iJtk At/Ax, for fij > , , 

- m,j, h + i n jVni ^ k At/Ax - Hjvm+t^kAt/Ax, for ^ < 1 ' 
j 



>k At + -a S) iVni im)k AiiAt - a/ va c<i n iJtk At r)^ k 

^li,m,j,k 



m=l 



m=l m=l 

where thj,* ~ n (%i, Hj,tk) is the number of neutrons at position X{ in direction Hj at time 
tk and a St i = <J s (xi). Also, r}fj mk and Vij,k are independent Gaussian 7V(0, l)-distributed 
numbers for each i, j, k, m. Notice that f fl3|) is an Euler-Maruyama approximation [13 1 120 1 121] 
to the system of Ito differential equations ([7]) and is a special case of the stochastic difference 
system (Q. 

The problem is solved numerically using two independent computational procedures, 
i.e, numerical solution of the SPDE is compared with Monte Carlo calculations. Equation 
( H~3|) is solved computationally with / = 80 equal intervals in position x and J = 40 equal 
angular intervals. The value chosen for the time interval is At = 0.125. In the Monte Carlo 
calculations, 1000 neutrons per second enter isotropically on the left side. Each neutron is 
followed individually in the Monte Carlo procedure with each neutron checked for a scatter, 
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a leakage, or an absorption at each time step of 0.1 seconds. Calculational results for 100 
sample paths using the two independent computational approaches are given in Table 2. 
The means and standard deviations of the calculated number of neutrons escaping from the 
left side and from the right side are given for the time interval t = 49 to t = 50. The two 
approaches agree well. In Figures |2] and |3j the calculated leakages for one sample path are 
compared for the two approaches. 



SPDE Calculation 



Monte Carlo 



800 



700 



800 




Figure 2: Calculated left leakages per second from time t = to t = 100 for one sample path 
using Monte Carlo and SPDE (flOl) . 



SPDE Calculation 



Monte Carlo 



1 50 



1 50 




1 OO 



Figure 3: Calculated right leakages per second from time t — to t — 100 for one sample 
path using Monte Carlo and SPDE ffTO]) . 



The second problem involves neutrons slowing down in a homogeneous medium. It is 
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assumed that at time t = 0.0, there are neutrons with energy between eV and 10 eV 
and 400 neutrons with energy between 10 eV and 20 eV. Also, there is present a constant 
source, Q(E), of neutrons where Q(E) = | q^^^q^ ^ < 10^ ^' 1US ' fo° Q(^) ^E = 220 
source neutrons are produced per second. For this problem, the source, Q(E), is assumed 
to be constant and does not vary randomly. It is furthermore assumed that the total 
cross section satisfies v(E)a(E) = 1.0 for < E < 20, the capture cross section has 

the form v(E M E) = ( ?» and the mattering cross section has the 

v ' v ' \ 1.0 for < E < 10 ta 

form v(E')a(E')f(E', E) = | o'o^for°0 <^E^ <W ^ "^° r convemence ' ^ ne P r °duct of the 
speed, v(E'), with the cross sections are given. The scattering kernel, a(E')f(E', E), is thus 
assumed to be piecewise continuous and proportional to l/y/W rather than, for example, 
to \jE' such as for a hydrogen-moderated system. Also, notice that the cross sections are 
consistent in the sense that 



/•20 

v(E')cr(E') = v{E')a c {E') + / v{E')a{E', E)f(E', E)dE=\ for all < E' < 20. 

Jo 

Furthermore, for this problem, it can be shown that the mean number of neutrons with 
energies between 10 eV and 20 eV is equal to 400 for time t > and the mean number 
of low-energy neutrons with energies between eV and 10 eV approaches 180 as time t 
increases. In this problem, the stochastic behavior of the number of neutrons with energies 
between eV and 10 eV is of interest. 

This problem is solved numerically using two independent computational procedures. 
Specifically, numerical solution of SPDE ( [121 is compared with Monte Carlo calculations. In 
the numerical solution of Equation (|12p . G = 20 energy groups of equal width AE = 20 /G 
are used. Equation (1121) is solved using the following stochastic difference equations at 
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discrete times tk = kAt where At = 0.02: 

G 

n gM+ i = n g>k + q g<k AEAt - v g cr g n g>k At + ^ v g :a g , f g > tg n g > >k AEAt (14) 

g'=l 

G 

-^Va~a~n~k~At r] { g c l - \ /v g a gfg,g' n g,k AEAt r]^ g } k 

g'=l 

G 

+ V v 9' a 9'f9',g n 9',k A EAt r$l )k 

g'=l 

where n g , k ~ n {E g , tk) is the number of neutrons in the gth energy group at time tk = kAt and 
a gfg,g' = a (E g )f(E g} E g :). Also, T)g k , rj g tr g ] k ~ A/"(0, 1) are independent normally distributed 
numbers with mean and variance 1 for each g,g',k. Difference system (1T4"]) is a special 
case of the stochastic difference system ([6]). In the Monte Carlo procedure, the neutron 
population in each energy group is checked at each time step for an absorption, an energy 
group change, or for an addition from the neutron source. This procedure is continued for 
each time step until the final time t = 2. 

Calculational results for 100 sample paths using the two independent computational 
approaches are given in Table 3. The means and standard deviations of the calculated 
number of neutrons with energies between eV and 10 eV and with energies between 10 
eV and 20 eV are given for time t = 2. The two different computational approaches agree 
well. In Fig. HI the calculated number of neutrons with energies between eV and 10 eV are 
given from time t = to t = 2 for one sample path for each calculational method. Again, 
the results are very similar for the two different calculational procedures. 

5 Applications 

Three possible applications of the stochastic neutron transport equations (jSj) are in develop- 
ment of new computational methods for solving neutron transport problems, in sensitivity or 
perturbation analysis, and in testing of computational and analytical methods. Numerical 
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Figure 4: Calculated number of neutrons with energy less than 10.0 from time t = to t = 2 
for one sample path using Monte Carlo and SPDE fTl2l . 

methodology for stochastic partial differential equations is advancing rapidly. It is probable 
that numerical techniques for solving stochastic neutron transport equations may eventually 
out-compete Monte Carlo techniques. Stochastic neutron transport equations can also be 
useful, for example, in perturbation studies to evaluate the effects of small changes in reac- 
tor parameters. In addition, solutions of stochastic neutron transport equations can provide 
independent checks on numerical or analytical approaches such as Monte Carlo methods. 

Consider a very simple example of applying stochastic neutron transport equations in a 
perturbation study. Consider energy- dependent transport in a homogeneous medium with 
the neutron density given by (fl2|) . Assume, for this example, that there is no neutron source 
and there are no collisions other than capture collisions where the capture cross section varies 
with time. The stochastic transport equation, for this problem, reduces to: 



dN(E,t) 



-va, 



(E,t)N(E,t) - y/va c (E,t)N(E,t) 



d 2 W(E,t) 



(15) 



dt v ' ' v ' ' v v ' ' v ' ' dEdt 

Suppose that a c (E,t) is perturbed to cr c (E,t) = a c (E,t) + Aa c (E,t) for t > 0. We wish 
to estimate, for the perturbation, the total number of neutrons as well as the change in the 
variability in this number with time t. In particular, if the perturbed neutron density is 
N(E,t) + AN(E,t), we wish to estimate the change in the total number of neutrons as a 
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function of time, i.e., An(t) = JJf° AN(E, t) dE. To derive equations for this quantity, let 
E k = kAE where AE is a small energy interval. Integrating ( fl5|) over the interval [Ef-, E k+ \], 
then 

— —— = -v k a c (E k ,t)n k (t) - \/v k a c {E k ,t)n k {t) — — — (16) 
for k — 0, 1, . . . where n k (t) = J E k+1 N(E, t) dE. In addition, using Ito's formula (e.g., [13]). 

= -2v k a c {E k ,t)n 2 k {t) + v k a c {E k ,t)n k {t) - 2n k (t)^v k a c (E k ,t)n k (t) ^|^- (17) 

From ( fl6|) and (fl7|) . expressions for E{n k {t)) and E{n 2 k {t)) are readily obtained and then, 
E(n(t)) = E(J2 k n k (t)) and £(n 2 (t)) = E((J2 k n k {t)) 2 ). Finally, expressions for E(n(t) + 
An(t)) — E(n(t)) and Var(n(i) + An(t)) — Var(n(t)) are derived as: 

£(n(t) + An{t)) - E(n(t)) = J N(E,0) (e~ ft _ e -j* Wc (£, s )<fej dE ^ 

and 

Var(n(t) + An(t)) - Var(n(t)) = 

- jT iV(£, 0)e" 2 ™^.')* (J va c {E, s)e% ^{E,r)dr ds \ dE ^ 

Therefore, using the stochastic neutron transport equation to analyze the effect of the pertur- 
bation for this transport process, not only are equations derived for estimating the average 
effect of the perturbation but equations are also obtained for estimating the change in the 
variability for the perturbation. Indeed, for this perturbation problem, equations (fTB|) and 
( 1191) clearly indicate that the mean change and the change in the variability are both pro- 
portional to the initial number of neutrons. 

6 Conclusions and Future Directions 



Stochastic difference and partial differential equations (SPDEs) are becoming increasingly 
important in applied mathematics [TTJ [151 EH 121] • In the present investigation, stochastic 
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difference equations and an SPDE are derived for neutron transport in a general three- 
dimensional medium. In the derivation procedure, the deterministic and stochastic terms 
in the differential equation system are simultaneously derived. First, a stochastic difference 
system is constructed. Next, an SDE system is derived. Finally, a particular SPDE follows 
from the SDE system. The stochastic difference equations and the SPDE for the neutron 
angular densities are given by © and dSJ), respectively. SPDEs for special cases of this equa- 
tion are given by ( TTOl) and (112p for one-dimensional plane geometry and for a homogeneous 
medium, respectively. The stochastic equations generalize the deterministic neutron trans- 
port equations and include random influences due to interactions and sources. Hence, certain 
random phenomena, such as fluctuations during reactor startup, can be studied using these 
SPDEs. The stochastic difference equations for neutron transport are solved numerically 
and compared with independently formulated Monte Carlo methods. The computational 
results between the two different numerical methods are in good agreement supporting the 
accuracy of the stochastic neutron transport derivation procedure. 

Future work may include appropriately extending the derivations of the present inves- 
tigation to include the random influence of prompt and delayed neutrons [9j [16j ITT] . In 
addition, stochastic difference equations and an SPDE may be developed to model the ran- 
dom behavior of neutron transport in spherical and cylindrical geometries. 
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TABLE 1 

Possible Changes in the Number n = n(xi, yj, z k , Hh (f> m , E g , t p ) for Time At 



Lyiiange ^AA/ij 


Description 


.riooaoiircy m lime laz 


li x v g n(xi-i, yj, z k ,ni,4> m ,E g , t p )At/ Ax 


T 

In 


a yz-tace [h x > 0) 


1 


-/j, x v g n(x i+1 ,yj, z k ,n h m , E g , t p )At/ Ax 


T 

In 


a yz-tace \\i x < 0) 


1 


II / J. n i \ A -t / A 

-\Vx\v g n{xi, yj, z k , ii h (pm, Eg, t p )At/Ax 
fj, y v g n(xi, yj-i, z k , Hi, 4> m , E g , t p )At/ Ay 




Out a yz-face 


1 


In 


a xz-face {Hy > 0) 


i 
i 


-fx y v g n(xi, y j+1 , z k , Hi, <f> m , E g , t p )At/Ay 


In 


a xz-face (Hy < 0) 


1 


-\Hy\v g n(xi, yj, z k , Hi, (fi m , E g , t p )At/Ay 




Out a xz-face 


1 


H z v g n{xi, yj, z k -i, Hi, 4> m , E g , t p )At/ Az 


In 


a xy-face (fx z > 0) 


1 


-li z v g n(xi, yj, z k +i, m, <j> m , E g , t p )At/ Az 


In 


a xy-face (Hz < 0) 


1 


-\Hz\v g n(xi, yj, z k , Hi, <f> m , E g , t p )At/Az 




Out a xy-face 


1 


-1 




Capture 


Pc 


-l/c 




Transfer out 


Ptrl 


1 




Transfer in 


Ptr2 


1 




Source 


PQ 
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TABLE 2 

Monte Carlo (MC) and SPDE Calculational Results for 100 Sample 
Paths for the Leakage for the Time Interval t = 49 to t = 50 



Average Number 
Out Left Side 


Standard 
Deviation 


Average Number 
Out Right Side 


Standard 
Deviation 


704.93 (MC) 
694.32 (SPDE) 


23.01 (MC) 
21.05 (SPDE) 


100.22 (MC) 
106.75 (SPDE) 


9.93 (MC) 
7.57 (SPDE) 
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TABLE 3 

Monte Carlo (MC) and SPDE Results for 100 Sample Paths for the Number of 
Neutrons at Time t = 2.0 With Energy Less Than 10.0 or Between 10.0 and 20.0 



Average Number 

With Energy 
Less Than 10.0 


Standard 
Deviation 


Average Number 
With Energy 
Between 10.0 and 20.0 


Standard 
Deviation 


156.83 (MC) 
156.97 (SPDE) 


11.44 (MC) 
10.28 (SPDE) 


399.97 (MC) 
400.92 (SPDE) 


14.39 (MC) 
13.42 (SPDE) 
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Figure Captions 



Fig. 1. A Brownian sheet on [0,5] x [0,5]. 

Fig. 2. Calculated left leakages per second from time t — to t — 100 for one sample path 
using Monte Carlo and SPDE ffTO]) . 

Fig. 3. Calculated right leakages per second from time t — to t — 100 for one sample path 
using Monte Carlo and SPDE (TT0]). 

Fig. 4. Calculated number of neutrons with energy less than 10.0 from time t = to t = 2 
for one sample path using Monte Carlo and SPDE ( fl2l) . 
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